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Summary 

Many in the commercial aviation industry consider 
wind shear to be a major safety concern. NASA has 
been developing airborne, forward-looking, system con- 
cepts to detect hazardous wind shear. Doppler radar 
and LIDAR are two candidate concepts being tested 
to provide this capability. An inherent limitation of 
this type of system is its inability to measure velocities 
perpendicular to the line-of-sight. Although these sys- 
tems can detect the presence of a wind shear by measur- 
ing the divergence of the horizontal wind profile, their 
inability to measure microburst downdrafts can result 
in a significant underestimate of the magnitude of the 
hazard. However, the vertical wind component can be 
inferred from the horizontal wind profile through the- 
oretical or empirical analytical models of microbursts. 

A simple analytical microburst model has been devel- 
oped for use in estimating vertical winds from horizon- 
tal wind measurements. It is an axisymmetric, steady 
state model that uses shaping functions to satisfy the 
mass continuity equation anti simulate boundary layer 
effects. The model is defined through four model vari- 
ables: the radius and altitude of the maximum horizon- 
tal wind, a shaping function variable, and a scale factor. 
The model closely agrees with a high-fidelity analytical 
model and measured data, particularly in the radial di- 
rection and at lower altitudes. At higher altitudes, the 
model tends to overestimate the wind magnitude, rela- 
tive to the measured data. 

Introduction 

Many in the commercial aviation industry consider 
wind shear to he a major safety concern. Wind shear 
is a rapid change in wind velocity (speed or direction) 
over a relatively short distance. It. can be found in 
a variety of weather conditions such as gust fronts, 
sea-breeze fronts, and mountain waves. Hazardous 
wind shear is most often associated with the convective 
outflows of thunderstorms. One such convective activity 
is the microburst, which is a particularly lethal form 
of wind shear when encountered at low altitude. 1 he 
microburst is a strong localized downdraft of 2 to 10 
minute duration, which causes a significant outflow as 
it impacts the ground. 

The hazard of wind shear arises principally from its 
deceptive nature. In a typical microburst encounter, 
the airplane initially experiences a performance increas- 
ing headwind as it penetrates the outflow region of the 
microburst. This will cause either an increase in the 
airplane’s airspeed or a deviation above the intended 
flight path, or both. Sensing this, the pilot may take 
corrective action to return to the intended flight path. 
'Phis is normally achieved by reducing thrust and pitch 
resulting in a reduction in the energy of the airplane. 


As the airplane enters the core of the microburst, there 
is a rapid decrease in headwind, shifting to a tailwind. 
'Phis is accompanied by a strong downdraft. The com- 
bination of these wind changes tends to reduce airspeed 
and push the airplane toward the ground. If the energy 
of the airplane is sufficiently low it may not be able to 
recover from the encounter. 

The Federal Aviation Administration, as lead 
agency for civil aviation safety, has established an in- 
tegrated wind shear program plan aimed at reducing 
the hazard of low-level wind shear. The program plan 
incorporates the expertise of industry, universities, and 
various government agencies such as the National Aero- 
nautics and Space Administration (NASA), the Na- 
tional Oceanic and Atmospheric Administration, and 
the Department of Defense into a multi-year research 
and development efTort. The objective of NASA s con- 
tribution to the wind shear program is to provide the 
technology base to reduce the risk of low-altitude wind 
shear through airborne detection, warning, and avoid- 
ance. Key elements of the NASA research effort in- 
clude characterization of the wind shear phenomena in 
the aviation context, airborne remote-sensor technol- 
ogy that provides forward-looking avoidance capability, 
and flight-management system concepts that provide 
timely and accurate transfer of weather hazard infor- 
mation to flight crews. The direction of the NASA re- 
search thrust is toward developing system concepts that 
embrace forward-looking sensor technology. Thus pro- 
viding the flight crew with awareness of the presence of 
wind shear with enough time to avoid the affected area 
or escape from the encounter. 

A fundamental requirement for a forward-looking, 
airborne wind shear detection system is the ability to 
estimate reliably the magnitude of the wind shear haz- 
ard that would be experienced by an airplane if it were 
to continue along the line-of-sight. The hazard is a 
function of both the horizontal and vertical wind com- 
ponents. Doppler radar and LIDAR are two candidate 
concepts being tested to provide this capability. Both 
systems measure the Doppler shift from aerosols, rain 
drops and other debris in the air to determine the line- 
of-sight relative velocity of the air. An inherent limi- 
tation of this type of system is its inability to measure 
velocities perpendicular to the line-of-sight. Although 
these sys tents can detect the presence of a microburst 
by measuring the divergence of the horizontal wind pro- 
file, their inability to measure the downdraft can result 
in a significant underestimate of the magnitude of the 
hazard. However, the vertical wind component can be 
inferred from the horizontal wind profile through the- 
oretical or empirical analytical models of microbursts. 
The simple analytical model described in this report 
was developed for this purpose. 
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Symbols 

r radial coordinate (distance from downhurst 
center), meters 

ti wind velocity in the x-direction, m/sec 

u T wind velocity in the r-dircction, in/sec 

u wind velocity in the y-dircction, m/sec 

w wind velocity in the z-direction, m/sec 

x longitudinal coordinate, meters 
y lateral coordinate, meters 

z vertical coordinate (positive up), meters 


Model Development 

The analytical microhurst model described in this 
report is a modification of the Osegucr a/ Howies 
downhurst model described in reference 1. The 
Oseguera/ Bowles model was developed as a simple ana- 
lytical microburst model for use in studying wind shear 
escape procedures. It is an ax [symmetric, steady state 
model that uses shaping functions to satisfy the mass 
continuity equation and simulate boundary layer effects. 
The following describes the derivation of this model and 
the subsequent modifications. 

The axisymmetric, steady state mass continuity 
equation expressed in cylindrical coordinates is: 


du T d w u r 

dr + dz r 
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This equation can be satisfied by solutions of the form: 


provided 


where 


u T = f(r)p(z) 
W = ;/(r 2 )</(;) 
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/(r) radial shaping function of the horizontal wind 
velocity, m/scc 

g(r 2 ) radial shaping function of the vertical wind 
velocity 

p(z) vertical shaping function of the horizontal wind 
velocity 

q(z) vertical shaping function of the vertical wind 
velocity, m/sec 


A scaling factor, 1/sec 


The shaping functions 
Oseguera/Bowles model were: 

selected for 

the 
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(8) 

q(z) = -A je - l] - 
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where 


R radius of downhurst shaft, meters 

z* characteristic height (out of boundary layer), 
meters 

e characteristic height (in boundary layer), 
meters 


Figures 1 and 2 show the characteristic shape of these 
functions. 

The resultant horizontal and vertical wind equations 
for the Oseguera/Bowels model were obtained by com- 
bining these shaping functions. 

Ur = [, _ e -0- 2 /« 2 >] - e-(*/ f )] (10) 

w = _Ae-( r2 / R2 ) {e [*"<*/*> - l] - z* [»-<*/*’> - l] } 

( 11 ) 

This model was initially considered for use in esti- 
mating the downdraft from measured horizontal winds. 
The downdraft is estimated by searching for the model 
variable values that yield the best fit of the model to 
the measured horizontal winds and then using these val- 
ues to compute the downdraft. The Oseguera/ Bowles 
model was found to be ill suited for this purpose. The 
radial shaping function of the horizontal wind f(r) was 
sufficiently different from the characteristic horizontal 
wind profile of a microburst that it was difficult to 
achieve a good fit to the measured winds. This then re- 
sulted in an inaccurate downdraft estimate. The modi- 
fication of the Oseguera/Bowles model involved defining 
a new /(r) that more accurately matched the character- 
istic horizontal wind profile of a microburst. The new 
shaping function was defined as: 
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where 

a shaping function variable 
f3 shaping function variable, meters 

The radial shaping function of the vertical wind 
g(r^) is derived from equation 4 which can be written 

, 2 , 2 awwi 
9<r > = 
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Figure 3 shows the characteristic shape of the new 
/(r) and </(r 2 ) functions. The vertical shaping functions 
p(z) and tj{z) were not changed. The new horizontal 
and vertical wind equations are: 
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Model Characteristics 

The radius of the peak horizontal wind ( r v ) at a 
given altitude (r) can be derived by setting the partial 
derivative of u witli respect to r equal to zero and 
solving for Tp. This results in an equation of the form. 



(&\ _i 
\(f* J 2 

(17) 

or 

r. 
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(18) 

Therefore 

0 = 2, \J 2r l n 

(19) 

The altitude of the maximum horizontal wind z m 
is dependent on the shaping functions p{z) and q(z). 
These functions were not changed from those of refer- 
ence 1, where it was shown that: 
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The velocities u r and w can be expressed in terms 
of r p 7 Zm . the shaping variable a, and the scale factor 
A by combining equations 19, 21 and 22 with equations 
15 and 16. 
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c 2 = -2.75 


(24) 

(25) 

(26) 


The equations for the wind along a radial can be 
normalized by dividing u r by the peak horizontal wind 
Up, and w by the downdraft at the center of the mi- 
croburst wo* 
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(30) 

From equations 27 and 28 it can be seen that the 
shape of the normalized form of the wind equations is 
only dependent on the value of the shaping variable 
o. Figures 4 and 5 show how a effects the shape of 
the normalized horizontal and vertical wind functions, 
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respectively. Increasing a tends to steepen the rate of 
decay of the horizontal winds and increase the updraft 
of the vertical winds (updraft = w/wq < 0). 

Figures 6 and 7 show how the normalized hori- 
zontal and vertical winds compare with other models 
and measured data. The TASS (Terminal Area Sim- 
ulation System) mean profile data is the average of 
eight axisymmetric microburst simulations taken at or 
near the time of maximum outflow velocity. TASS is 
a numerical, time-dependent, multi-dimensional, non- 
hydrostatic cloud model which has been used exten- 
sively in the study of microbursts (refs. 2 and 3). The 
'PASS vertical wind data presented in figure 7 was taken 
at an altitude of 500 meters. The JAWS (Joint Airport 
Weather Studies) mean profile shown in figure 6 is radar 
measured data taken from reference 4. 

The Oseguera/Bowles model, shown in figures 6 
and 7, considerably overestimates the horizontal winds 
while underestimating the vertical winds. This was 
the primary reason for developing the new shaping 
function. With a = 2, the new model more closely 
models the horizontal wind decay while avoiding the 
excessive updrafts which result from larger values of <\. 

The altitude dependent shaping functions, p(z) and 
q(z), were not changed from those used in the Os- 
eguera/Bowles model. Therefore the vertical profiles of 
the normalized horizontal and vertical winds were un- 
changed. Figures 8 and 9 sltow a comparison of these 
profiles with radar measured values. Presented in figure 
8 is the peak horizontal wind at each altitude (u ?> ) nor- 
malized by the overall maximum horizontal wind (w m ), 
plotted against the altitude (z) normalized by the alti- 
tude of maximum horizontal wind (z m ). The equation 
for the model derived u m is: 

= -c r »)e< , /2«> (31) 

Presented in figure 9 is the vertical profile of the 
downdraft at the core of the microburst (u>o) normalized 
by the core downdraft value at the altitude of maximum 
horizontal velocity (?ey m ). The equation for the model 
derived u>om is: 

wo™ = -A (— [c Cl - 1] - ^ [c* - 1]} «<*/“> (32) 

I Cl C 2 ) 


tildes. At higher altitudes, the model tends to overes- 
timate the wind magnitude, relative to the NIMROD 
data. 

This model was developed to estimate the down- 
draft from horizontal velocity measurements. However, 
it could also be used in a similar manner as the Os- 
eguera/Bowles model, as a simple microburst simula- 
tion model. The equations required to use this model 
for microburst simulation are provided in the appendix. 

Concluding Remarks 

A simple analytical microburst model has been de- 
veloped for use in estimating vertical winds from hor- 
izontal wind measurements. It is an axisymmetric, 
steady state mode] that uses shaping functions to satisfy 
the mass continuity equation and simulate boundary 
layer effects. The model is defined through four model 
variables: the radius and altitude of the maximum hor- 
izontal wind, a shaping function variable, and a scale 
factor. The model closely agrees with a high-fidelity 
analytical model (TASS) and measured data (JAWS, 
NIMROD), particularly in the radial direction and at 
lower altitudes. At higher altitudes, the model tends to 
overestimate the wind magnitude, relative to the NIM- 
ROD data. 

References 

1. Oseguera, Rosa M.; Bowles, Roland L.: A Simple, 

Analytic 3- Dimension al Downburst Model Baaed on 
Boundary Layer Stagnation Flow. NASA TM- 100632, 
July 1988. 

2. Proctor, F. H.: The Terminal Area Simulation System, 
Volume I: Theoretical Formulation. NASA CR-4046, 
April 1987. 

3. Proctor, F. II.: The Terminal Area Simulation Sys- 

tem, Volume !l: Verification Cases . NASA CR-4047, 
April 1987. 

4. lljelmfelt, M. R.: Structure and Life Cycle of Microburst 
Outflows Observed in Colorado. Journal of Applied 
Met ecology , Vol. 27, 1988, pp. 90. 

5. Fujita, T. Theodore: Tornadoes and Downbursts in the 
Context of Generalized Planetary Scales. Journal of the 
Atmospheric Sciences , Vol. 38, August 1981, pp. 1511- 
1534. 


The NIMROD (Northern Illinois Meteorological Re- 
search On Dow id Hirsts) data presented in figures 8 and 9 
are Doppler radar measured winds of a microburst event 
presented in reference 5 (figure 30). The microburst 
translational speed of 15 m/s was subtracted from the 
horizontal wind speed. 

The model agrees closely with the measured data 
presented in figures 8 and 9, particularly at lower alti- 
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Appendix 


Velocity Equations 


Micfoburst Model Simulation Equations 

Provided below arc equations for the wind velocities 
and their spatial derivatives at the point P in the 
coordinate system shown in sketch a. The microburst 
center is located at the origin of the axis system. 



Sketch a. Coordinate syslern. 


Model Variables 


This model can be fully defined by specifying the 
radius, altitude and magnitude of the maximum outflow 
( Vp , z 77i and u m , respectively), and the value of the 
shaping variable a (o = 2 is recommended). The scale 
factor A can then be determined from: 
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Spatial Derivatives 


— = - L=l(*/=m) 

flx 2 


x 2 (x 2 + y 2 ) 

.. 2 or 


Or — 1 1 




x e i 


(A 7) 


Ou -Xxy (j 2 +y 2 ) a 1 fc,(.-/ 2m )_ e cj( 2 / Jm ) 

<hf 2r p 2 « I 


/>« 



(2I8) 


x e 


2-(*W)°/r p 2 


(219) 


5 



dw A z(x*+y*) a ~' \ , , (* 2 + yT ~ 

dx r p 2ct ^ 2 r p 2 " 
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Increasing Altitude 
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Increasing Radius >■ 

Figure 1 . Characteristic variation of radial shaping functions. 



Figure 2. Characteristic variation of vertical shaping functions. 
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